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Abstract — General formulas of entanglement concentration are 
derived by using an information-spectrum approach for the 
i.i.d. sequences and the general sequences of partially entangled 
pure states. That is, we derive general relations between the 
performance of the entanglement concentration and the eigen- 
values of the partially traced state. The achievable rates with 
constant constraints and those with exponential constraints can 
be calculated from these formulas. 

Index Terms — Information spectrum, Entanglement concen- 
tration, Exponents, Maximally entangled state 

I. Introduction 

VARIOUS quantum information processings are proposed, 
many of which require maximally entangled states as 
resources, e.g., quantum teleportation and dense coding etc[2], 
[1], [3]. Hence, it is often desired to generate maximally 
entangled states. However, the realized state is not necessarily 
a maximally entangled state. Thus, entanglement concentration 
is used for producing maximally entangled states (MES) from 
partially entangled pure states only by local operation and clas- 
sical communication (LOCC), while entanglement distillation 
is used for producing them from partially entangled mixed 
states by LOCC. Therefore, entanglement concentration is an 
important issue in the field of quantum information. 

In information theory, we often assume that the system is 
prepared as the independent and identical multiple copies of 
the given state. Such a condition is called independently and 
identically distributed (i.i.d.) condition. Under this condition, 
Bennett et al.[6] showed that the amount of entanglement of 
a partially entangled pure state |$)($| is described by the 
entropy H(p) of its partial traced state p :— Tr-^ B |<J>)(<1>|, 
which is called the reduced density matrix. That is, they proved 
that an MES with size 2 nH ^ can be asymptotically produced 
from n identical copies of the state |<&) (<J>| with a high enough 
probability. Furthermore, independently of the form of |<&) (<J>|, 
Hayashi and Matsumoto constructed a protocol satisfying the 
above property, which is called universal [7]. 

However, in the correlated physical system, the state of the 
total system cannot be regarded as independent and identical 
copies of a given state. In such a case, we have to treat 
general partial entangled pure state between two distinct 
parties. Indeed, this model is not so unnatural because the 
state on the total system is pure when this system is isolated 
from the other system. In this paper, as a general asymptotic 
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method to treat this model asymptotically, we focus on the 
information spectrum method and apply it to entanglement 
concentration. The information spectrum method has been de- 
veloped by Han and Verdu [9] for discussing general sequence 
of information sources/channels, and been established as a 
unified method to information theory in Han's textbook[13]. 
Indeed, this method has been applied to quantum information 
theory, for example, to quantum hypothesis testing[14] and 
quantum channel coding [15]. In this paper, we apply this 
method to entanglement concentration, and characterize the 
asymptotic production rate of a general sequence of partially 
entangled pure states without any assumption. The information 
spectrum method used in this paper is slightly different from 
the original Han-Verdu's method, and is close to Nagaoka- 
Hayashi's method[14]. 

In the derivation of our general asymptotic formulas, 
we essentially use the majorization method established by 
Nielsen[4]. Based on this method, he developed a necessary 
and sufficient condition for the possibility of transforming 
from a partially entangled pure state |<J>i)(<J>i| to another 
entangled pure state |$2)($2| by using LOCC between the 
two parties Ha and Hb- This condition is characterized 
only by the eigenvalues of their reduced densities pi := 
Tr HB |$i)<$i|, (» = 1,2). 

Moreover, even in the i.i.d. case, the knowledge of the 
asymptotic production rate is not sufficient for estimating the 
production rate of MES for a given finite number of copies. In 
channel coding or source coding, for this analysis, we usually 
focus on the error exponents, i.e., the exponential rate of error 
probability because the error goes to exponentially when 
we choose our code suitably. In entanglement concentration, 
when we fix the production rate to a constant number less 
than the entropy rate, the optimal failure probability goes to 
exponentially. Hence, based on its exponential rate (failure 
exponent), we can roughly estimate the failure probability 
for a given finite number of copies. As preceding researches, 
Hayashi et aZ.[8] derived the failure exponent of entanglement 
concentration in the i.i.d. case based on the method of types, 
and Hayashi and Matsumoto [7] did that of their universal en- 
tanglement concentration protocol. In this paper, we calculate 
the failure exponent of entanglement concentration in a more 
general setting. 

In most problems in information theory, in the i.i.d. case, 
the correct (or success) probability exponentially goes to 
when the rate is strictly better than the optimal rate. This 
exponential rate is called the correct (or success) exponent, 
and is one of famouse issues in information theory. Hayashi 
et aL[8] and Hayashi and Matsumoto [7] treated the success 
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exponent of entanglement concentration in the i.i.d. case. This 
paper proceed to the general sequence of partially entangled 
pure states. 

One may think that such an exponential treatment is not 
essential. It is, however, more difficult to obtain the error 
and correct exponents asymptotically and tightly than the 
asymptotical optimal production rate. Hence, in order to derive 
these tight bounds of both exponents, we need better and 
more simple non-asymptotic evaluations. That is, such a non- 
asymptotic evaluation should be a better and more simple 
approximation for the optimal value. Therefore, even though 
the optimal correct exponent is useless, the non-asymptotic 
evaluations used for its derivation is quite useful. 

Furthermore, the optimal rates with exponential constraint 
are characterized by Renyi entropy in the i.i.d. case. In this 
paper, we derive the same formulas under a weak assumption 
for the Renyi entropy. Using these formulas, we characterize 
the optimal rates based on the partition function. 

Finally, we have to explain our formulation of entanglement 
concentration. There are two formulations in source coding. 
One is fixed length, in which the coding length is fixed, i.e., 
is independent of the input data. The other is variable-length, 
in which the coding length is variable, i.e., depends on the 
input data. Similarly to source coding, we can consider two 
similar formulations in entanglement concentration. In Bennett 
et al. [6]'s protocol and Hayashi and Matsumoto[7]'s protocol, 
a local measurement is required as the first step, and the length 
of the MSE generated finally depends on the data of this 
local measurement. Hence, their protocol is a variable-length 
entanglement concentration. 

On the other hand, based on Nielsen's result[4], Hayashi 
et al. [8] discussed entanglement concentration protocols pro- 
ducing the MES with the fixed size. Hence, such protocols 
are called fixed-length entanglement concentration, which are 
classified into two formulations as follows. In the first for- 
mulation, we produce, without a failure, an approximately 
MES from a partially entangled pure state. Its performance is 
represented by the size of the MES and the fidelity between the 
appropriate MES and the final state. This kind of entanglement 
concentration is called deterministic fixed-length entanglement 
concentration (DFLEC). In the other formulation, we produce 
an MES itself, allowing a failure probability, from a partially 
entangled pure state. The performance of this protocol is eval- 
uated by the size of the MES and the failure probability. This 
protocol is called a probabilistic fixed-length entanglement 
concentration (PFLEC). Hayashi et a/. [8] treated these two 
formulations in the i.i.d. case. In this paper, we discuss them 
in a more general model. 

This paper is organized as follows. In section|ll] we give the 
mathematical definitions of the optimal rates with respective 
conditions, (constant constraint, exponetial constraint) for the 
genereal sequence of partially entangled pure state in two 
formulations of FLEC. As the main results, characterizations 
of these quantities are given based on information spectrum. 
That is, we discover a general relation between the perfor- 
mance of entanglement concentration and the eigenvalues of 
the reduced density (partially traced state). In section |lll] the 
optimal rates of FLECs are characterized by the Renyi entropy. 



In section IIVI we apply these formulas to the case when the 
reduced density is given as a thermal state. In section [V] the 
performances of the two FLEC types in a non-asymptotic case 
are characterized by applications of Nielsen's result [4] and 
Lo and Popescu's results[17]. In section fVll the main result 
is verified by applying several lemmas described in section [V] 
to an asymptotic case. In section IVHI the relation between 
entanglement concentration and random number generation 
is discussed. The appendix [X] summarizes relations for the 
quantum analogue of the information spectrums based on the 
original definition [14], which are necessary for verifying the 
main result. 

II. Main results 

When the two distinct parties, Alice and Bob, have their 
respective systems Ha and Hb, the total system is described 
by the tensor product space Ha^Hb- In quantum information, 
as is mentioned in section [I] one of main issues is the char- 
acterization of entanglement between these distinct parties. If 
the state on total syste is a pure state $ 6 Ha <8> Hb, it 
is known that its entanglement between two parties can be 
characterized by the reduced density (partially traced state) 
p := Ttb |$) ($|. In particular, if the reduced density p is the 
completely mixed state -j^I, it is called maximally entagled, 
where cIa denotes the dimension of the system Ha- Hence, if 
the pure state ^ £ Ha <8> is maximally entangled, there 
exist completely orthogonal basis {e^} and {e^} on Ha and 
Hb, respectively such that 

While any quantum operation is mathematically described by 
trace-preserving completely positive (TP-CP) map, in the en- 
tanglement concentration of the initial pure state $ on distinct 
two parties Ha and Hb, our operation is often restricted to 
a quantum operation with an LOCC implementation between 
Ha and Hb - Hence, a deterministic fixed-length entanglement 
concentration (DFLEC) is an LOCC quantum operation C 
together with a maximally entangled state on a subspace 
H' A ®H' B , i.e., it is described as (C, ^). Since this protocol 
(C, <I>) transforms the initial pure state |$)($| to the final 
state C($) := C(|<1>) (<1>|), its performance is evaluated by the 
fidelity {^\C($)\%) and the size L(*) of which equals 
H(Tr B |*)<¥|). 

For a rigid analysis of the probabilistic fixed-length en- 
tanglement concentration, we have to discuss a measuring 
operation that describes a quantum measurement with the final 
state as well as the probability distribution of the measured 
data. The measuring operation is given as a CP map valued 
measure / = {Ii}i whose sum is a TP-CP map; i.e., every 
is a CP map, and J^. /j is a TP-CP map. It is often called 
an instrument. When we perform a quantum measurement 
corresponding to / = on the system with a state p, we 
obtain the measured data i and the final state ^jrr^Ii(p) with 
the probability Tr/j(p). Hence, a probabilistic fixed-length 
entanglement concentration (PFLEC) of an initial pure state 
$ £ Ha 8> Hb is a two- valued instruments / = {Iq, h} with 
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an LOCC implementation satisfying that Zi($)/Tr/i($) is a 
maximally entangled state |^')( v l'| on a subspace Ti' A ® H' B , 
where is abbreviated to That is, the event 1 

corresponds to success, and the event does to failure. Thus, 
its performance is characterized by the failure probability 
Tr J ($) and the size L(I) :— of the final maximally 

entangled state. 

Here, we briefly discuss the relation between two kinds of 
fixed-length entanglement concentrations. For any PFLEC I = 
{Io,h} of $, the pair (h + 1 , Ji($)/ Trli($)) becomes a 
DFLEC and its fidelity between the final state and the desired 
maximally entangled state 7i($)/ Tr7i(<I>) is greater than the 
success probability of the DFLEC I = {Iq,Ii}: 



Tr 



+*>)(*) 



Tr/i($) 



> Tr/i($). 



(1) 



That is, for any a given PFLEC protocol, there exists a DFLEC 
protocol whose performance is better than the given PFLEC 
protocol. 

In the quantum system, if n systems are prepared identically 
to the system Ha <£> Hb, the total system is described by 
Hf 1 <g> H% n . If the state of every system H A <£> H B is the 
pure state $ and if each system is independently prepared, 
the state of the total system is written by the tensor product 
pure state Such a case is called the i.i.d. case. However, 
even if the state of each system coincides with each other, 
if they are not independent of each other, the state of total 
system is not a tensor product state. In order to treat such a 
general case, we focus on a general sequence of the pair of 
the joint system with distinct two parties Ha.u and Hb,u and 
the partially entangled pure state $„ e Ha,u <8> H-B,n with 
an asymptotic situation. Note that, in this notation, the space 
H.A,n and Hb,u are generalizations of 7i®" and H B m , and 
<!>„ is a generalization of the n-tensor product vector (J® 71 . 

In order to discuss the asymptotic optimal performance in 
such a general case, we optimize the production rate of MES 
with three asymptotic constraints for the failure probability or 
fidelity. Concerning the PFLEC, we focus on the following 
conditions: 

• Constant constraint: The asymptotic failure probability is 
less than a fixed constant. 

• Exponential constraint for the failure probability: When 
we choose a good DFLEC protocol, failure probability 
goes to exponentially. Hence, as another criterion, we 
restrict our DFLEC satisfying that the exponent of failure 
probability is greater than a fixed exponent. 

• Exponential constraint for the success probability: If we 
choose a bad DFLEC, the success probability goes to 
exponentially. Among such PFLEC protocols, if this 
exponent, i.e., the success exponent, is greater, the proto- 
col is worse. Hence, we can consider the optimization of 
the production rate of MES with the constraint that the 
success exponent is less than a fixed exponent. 



Thus, concerning PFLEC, we focus on the following values: 

log£(I") 



Bp(e) := sup < lim 
{J „ } L 



B e p(t) '■= sup \ lim 
{J „ } 1 



n 

log L{I n ) 



B* P (r) := sup < lim 



. log L(I n ) 



!Si7 "($ n ) < e} 
lim — log Tr >rj 
<r}. 



lim log Tr I 

n 



In the DFLEC case, we obtain several criteria by replacing 
the success probability in the above discussion by the fidelity. 
That is, we can define the following values: 

B D (e):= sup (lim - log L($ n ) 
{(c-,*„)} L n 

lim(*„ | C n ($„)|¥„) >l-e} 
sup (lim- logL(*„) 

(C",*„)} L n 

lim — log(l- (*„|C7 n ($ n )|*„>) >rj 
n j 

sup (lim- logi(* n ) 

Hm — log(* n |C"($ n )|* n ) <r). 
Hence, it is trivial from ([0 that 

B 1 (e)>B 2 (e), B e , D (r) > B e>2 (r), B* e<D (r) > B* e>2 {r). 

(2) 

In this paper, we treat a quantum analogue of information 
spectrums to analyze the above values. For such an analysis, 
we need the following definitions. For a self-adjoint operator 
X, we can denote the projection ^ > c Ei by {X > c}, 
where the spectral decomposition is given by X = ^\ XiEi. 
We can define the projections {X > c},{X < C}, {X < 
c}, etc. in a similar manner. Let p n be the reduced density 
Tr WB,„ and define 

K{a) := IurlTrp„{p„ - e~ na > 0} 

( c (a) :=ljm — logTrp n {p n -e- na > 0}. 
— n 

When the limit 



lim — \ogTi p n {p n - e~ na < 0} 
n 



(3) 



exists, we denote it by These definitions can also be 

written as 



(4) 



K(a) 










= limpn < 


— logPn.i 


< 






I n 




a J 


C(a) 










= lim — 


logPn 1 — 


log 


;p 




n 


I n 






C(«) 










= lim 


\ogp n \ 


log 


,p 


n 









(5) 



(6) 



where every p„ , is an eigenvalue of p n and can be regarded 
as a probability distribution. Hence, the quantity K(a), C°( a )> 
and C(a) denotes the degree of concentration of the e na - 
dimensional subspace. Note that the function C c ( a ) decreases 



4 



monotonically, while the function £(a) increases monotoni- 
cally. Indeed, in the classical case, the value K(a) gives the 
asymptotic performances of fixed-length source coding [12] 
and uniform random number generation [19], [13] with asymp- 
totic constant constraint. Moreover, the quantities C c ( a ) an d 
£(a) gives the asymptotic optimal performance of source cod- 
ing with the exponential constraint[12] and that of simulation 
of random process with KL divergence criterion [18]. As is 
mentioned in section rVTTI £(a) gives the asymptotic optimal 
performance of intrinsic randomness with KL divergence 
criterion [22]. 

As is mentioned in the following main theorem, the optimal 
production rate of MES can be characterized by how densely 
the eigen values of the reduced density matrix concentrate a 
small space. 

Theorem 1: Without any assumption, for every e € [0, 1] 
we have 

B D {e) = B P {() = sup{R\K(R) < e} 

R 

B e . D {r) = B e , P (r) = sup{i?|C c (i?) > r}. 

R 

WIuti the limit Q exists and there exists a real number a such 
that ((a) < ( c {a), we have 



B e,D( r ) =sup^a-r 



inf C(a')-- 



a' <a \ + -<r 



sup <! | -inf U(a')-j 



a < a 



a'<a} + -<r 



s; P (r)=sup{a-C(a)|C(a)<r}. 

a 

This theorem is proved in section IVII after preparing the 
appropriate discussion. 

Remark 1: As is mentioned in Nagaoka and Hayashi[14] 
the quantum versions of K(a),( c (a), and ((a) give the 
asymptotic performances of fixed-length source coding. In 
particular, the optimal rate with the constraint for the constant 
error exponent is given as 



sup{a- ((a)\({a) < r}, 



(7) 



which is almost similar to B* p (r). For a proof only of the 
classical case, see Han [12]. For a proof in the classical and 
quantum case, see Nagaoka and Hayashi [14]. 



III. Asymptotic formulas based on Renyi entropy 

In the classical and quantum fixed-length source coding of 
i.i.d. information source, it is known that the optimal rate with 
the constant constraint for error exponent is described by the 
Renyi entropy i/j(s) := log pf [21]. Concerning FLEC of 
the i.i.d. source, as is described in Theorem |2] Hayashi et 
al. [8] showed that this kinds of optimal rates can be described 
by the Renyi entropy. In this section, using Theorem Q we 
derive the same formula in a more general setting. 



Theorem 2: Hayashi et al.[8] When p n — p® n , the rela- 
tions 

B D {e) = B P {e) = H(p), Ve such that 1 > e > (8) 
B e , D (r)=B e , P (r)=sup r ^ {s) 

s>l 1 - 8 



KAr) = mm 



sr + ift(s) 



0<s<l 1 — S 

sr + i/j(s) 

B *e,D{r) = { 0<"<1 1 - S 

2^(i)+r 



(9) 
(10) 

if><-±V'(§)-^(|) 
otherwise 



(11) 



hold, where 

H{p) := -Trplogp, ^(s) :=logTrp s . 

In particular, the above formulas of some special cases are 
written as 

B e . D (r) = B e . P (r) =H 00 ifr>H 00 = lim -^'(s) 

s — >oo 

B* P {r) = V(0) if r > -V'(0) - ^(0), 

where 



lim 



s^oo s 

The following is the generalization of the above theorem. 

Theorem 3: Letting ip n (s) :— logTr p^, we assume that the 
limit -ip(s) := lim„ ^fcalsl exists and that its first and second 
derivatives ip (s) and ip (s) exist for s e (0, 1)U(1, oo). Then, 



H- <B D (e)=B P (e)<H+ 
B e , D (r) = B e , P (r) = sup ^f^ 

S >1 1 - S 

B eP (r)= mm — 

' 0<s<l 1 — s 



(12) 
(13) 

(14) 



2^(1) + 7 



ifr<-p'(I)-^(I) 



otherwise, 



(15) 



where 

H- := -^'(1 + 0), H+ := -^'(1 - 0). 
In particular, we have 

B e . D {r) = B e , P {r) = H 00 xfr>H oc = lim -^'(s) 
B* eP (r) = ^(0) if r > -^(+0) - ^(0), 



where H r 



lim., 



The equations (|9j, dlOi . and il It follow from the equa- 
tions Jl 31 . (I14> . and d!5l >. The equation (jSJl follows from 
the equation (I12l >. Hence, Theorem [3] can be regarded as a 
generalization of Theorem |2 Since Hayashi et al. [8]used 
the method of type, they proved Theorem |2] only in the 
finite-dimensional case. Hence, its infinite-dimensional case 
is proved by this paper first time. 
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Note that 



Remark 2: Under the same assumption as Theorem^ we From A25I . the equation /i(s r )' = holds. The derivative of 

can similarly prove that the numerator of /{ (s) is 

sup{<z-C(a)|C(a) <r}= min sr + ^ s \ - *) + r + ?(s)Y = ^'(s)(l - a) < 0, 

q Wl W 0<s<l 1 — S V ' 

which gives the optimal rate with the constant constraint for the final inequality inequality follows from G3. Therefore, 

error exponent in the fixed-length source coding. B e ^ D (r) = B e ^ P (r) = fi(s r ) = max s >i /i(s). 

Proof: As is discussed in Appendix E Gartner-Ellis Next ' we P™ e G3 for the case in which r > H x . From 

theorem [20] yields that EJ- if a > H °°> then C{a) _< r. Otherwise, C c (a) > r. 

Thus, B e £>(r) = B e p(r) = Hoo. Since the numerator of 

_ if a_<H + /J ( a ) equals 

sup (i _ s)a _ ^ (s ) > o if H+ < a < -iP (0) . 

o<s<i r + ip (s)(l - s) + ip(s) = r- ( c (-ip (s)) > 0, 

(16) we obtain f[(s) > 0. Therefore, 

C c (a\ = i supfl — s)a — ip(s) > if Hoc < a < H_ sup = lim = H^. 

2- 1 I 1 S >1 S >1 1 - S s^oo 1 - S 

oo if a < Hoc. — i 

ryj\ Proceeding to (I14t for the case in which r < — V' (+0) — 

ip(0), we define a r and s r by £(a r ) = F and s r :— s(a r ). 

Thus, we have 

CHft+o) - o) = -?'(+o) - ?(o) (is) (1 _ Sr)flr _ ?(sr) = r (26) 

^+0)=^. (19) -(1-^)-^)=, (27) 

Moreover, it follows from the discussion in Appendix 151 that __ 

—ft \ ■ c - -r» t \ ■ t f r /„ ,,,,/, \ Using ypj, we can calculate B p [r): 

ip(s) is convex. Since tp (s) exists tor s S (0, 1) U (1, oo), we e, - r 

have p* , v s r r + 4>{s r ) 

—ii B e P\ r ) =a r -r = . 

ip (s) > se(0,l)U(l,oo). (20) ' i-Sr 

For any real number a satisfying < a < -^(+0), we The derivative of the function / 2 (s) := £I±fcl(0 < s < 1) 

define s(a) by is g iven b Y 

a = -^'(s(a)). (21) ,// \ _ ^' (s){l - s) + r + ^(s) 

h[ s ) - (1-s) 2 

Hence, equations H61 and d!7l > yield that 

— From (1271 . the equation / 2 (s r )' = holds. The derivative of 

fin} = J _ a — + _, the numerator of fo(s) is given by 
tU I (l-«(o))o-^(*(o)) if < a < -V (0) 

(22 ) + =/(s)(l-s)>0 

f ifi?_<a because of @DJ. Therefore, S* p (r) = / 2 (s r ) = 

C c (a) = < (1 — s(a))a — ip(s(a)) if H QO <a<H- min s > 1 / 2 (s)- 

t oo if a < Hoc. Next, we prove i ll-It for the case in which r > —tp (+0) — 

( 23 ) ^(0). If a < -if(+0), then ((a) < r. Otherwise, 
First, we prove O for the case in which r < H^. In this ^ > ^_T hus - " follows from d that ^(r) = 

case, we can define a r and s r by ( c {a r ) = r and s r := s(a r ). lim 1 _ ++o (-i/'_(+0) - e) - C(-V> (+0) - e) = (+0) - 

Thus, we have (—i/j (+0) — tp(0)) = ip(0). Since the numerator of / 2 (s) is 



(1 - s r )a r - ip(s r ) = r (24) r + $ (s)(l - s) + ^(s) = r - C C (-^ (*)) > 0, 

-(1 - s r )i/j (s r ) - ip(s r ) = r. (25) then /^(s) > 0. Therefore, 



Using J24i . we can calculate B e jj(r) and B e p(r) as sr _|_ ^/,c s ~) sr _|_ ^c s ~) _ 

_' min — — lim — = ib(0). 

r + ^( Sr ) o<s<i 1-s s^o 1-s 

£>e d(t) — B E p(r) — a r = . __ 

1 — s r Next, we prove d!5i . We can calculate the derivative of 

The derivative of the function /i(s) :— life) (s > l) is given ^ 2 
by 



(C(a) - |) = 1 - s(a) - s'(a)a - ip' '(s(a))s'(a) - ^ 



ft/ \_ tp (s)(l - s) + r + ip(s) 1 i 

h( s ) — ^ _ g N 2 ■ =1 - s(a) - s'(a)a + s (a)a - - — - - s(a). 



6 



This derivative is if and only if s(a) - 
The second derivative is calculated as 



r, i.e., a 



> 0, (28) 



where the final equation follows from 1 = ip (s(a))s'(a) 
which can be derived from i2U . Thus, the function a i— > 
((a) — % is strictly convex, and its minimum value equals 



c(-^ (*)) + #(*). 

Hence, we have 



which is attained at a = —tp (i^ 



1-^(1) ifa>-V(i) 



where we use the equation c( — ^ (i)) 
which follows from (|22}. Since C(-?'(l/2)) = (§) - 
5^ (5), we have 

sup {a — r\((a) < r} 

a<-J/(l/2) 

\a r -r if r <_^(|)-p'(|) 
\ -^'(1/2) -r if r > -</> (i) - (|) ■ 

Remember that a r is defined such that C( a r) = r. Moreover, 
we have 



sup < a — r 

a>-^'(l/2) 



0_ if r <_^(i)_i^(i) 
2V>(±)+r if r >-^(|)-^'(i). 



Therefore, 



s e, d( v ) =sup|a-r 



inf C(a') 



a < a f- + — < r 



max< sup {a ~ r \C(a) < r } , 

[a<-^'(l/2) 




sup ^ a — r 

a>-^'(l/2) 

= fa L -r ifr<-^(|)- (±) 

\2^(|)+r if r > -i/> . 

Using a discussin similar to J14-I - we can show H5\ . ■ 

IV. Correlated system 

In this section, we consider the application of Theorem|3]to 
correlated systems. As an example, the initial state is assumed 
to be a ground state with the Hamiltonian J2i Hi + Hi^+i 
on the system (Ha <8> Hb)® 11 , where Hi is the Hamiltonian 
of the i-th joint system between A and B, and is 
its interaction term between the z-th and i + 1-th systems. 
However, it is not so easy to calculate i/j(s) in this case. Hence, 
we focus on a more ideal case. 

Assume that the total system (Ha <8> Hb)® 11 is isolated 
from other systems. We also assume that the system H!^ n 
is sufficiently large, and the interaction between the system 



Hf 1 and the system H% n is ideal so that the system H% n 
can be regarded as the heat bath of the system Tt® n . Now, 
we suppose that the Hamiltonian ^\ Ha^ + on the 

system 7i? n . Hence, the state of the total system is pure, 
and the reduced density on A is the thermal state with the 
Hamiltonian J2i Ha^+Ha^j+i. Now, we define the partition 
function as 

3(/3) :=lwi~\ogTrexp((3yHA,i + H A , iti+1 ). (29) 



Thus, when the inverse temperature is /3q and the partition 
function is continuous and differentiable, the tjj(s) can be 
calculated as 

4>{s) = hm - log Tr — - 

n VTrexpf/jo^-ffAi + ^u+i) 

= lim- logTrexp(s/3 o y^i?A,i + H A ,i.i+i) - s3(/3 ) 
n ' 

i 

=E(s(3 )- s E(p ). 
Hence, 

B D {e) = B P (e) = -po&(Po) + 3(/3 ) 



B e ,D{r) = B e ^ P (r) = sup 



r + E(s/3 ) - sE(p ) 



B* P (r) = min 

' 0<s<l 



mm 



s>l 1 - S 

sr + E(s(3 )- S E{f3o) 
1-s 

. sr + 3(s/3 )-sS(/? ) 



1 - s 



if r < ri/2 
23 - 3 (/?o) + r otherwise, 



where 



Note that the above formulas are based only on the partition 
function. Hence, it is expected to apply them to other cases. 
Moreover, we can derive similar formulas concerning classical 
and quantum fixed-length source coding. 

V. Non-asymptotic theory 

In order to derive general asymptotic formulas based on 
the quantum information spectrums, we need to prepare ap- 
proximate formulas for non-asymptotic setting based on the 
form of the reduced density p. For this purpose, we focus 
on majorization, because it gives a necessary and sufficient 
condition for the possibility of transforming from a partially 
entangled pure state to another entangled pure state 

|$2)( < i ) 2| by using LOCC between the two parties Ha and 
W B [4]. Suppose that p = (p x , . . . ,p d ) and q = (q x , ...,q d ) 
are probability distributions. The probability p majorizes q, 
(equivalently q is majorized by p), written p y q, if for each 
k in the range 



E*i- 



The elements indicated by J, are taken in descending order; 
for example, p[ is the largest element in (pi, ■ ■ ■ ,Pd)- The 
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majorization relation is a partial order. To discuss entangle- 
ment transformation, we need to treat probability distributions 
consisting of eigenvalues of a reduced density p. The reduced 
density p majorizes another reduced density a written p y a, 
if the probability distribution p(p) consisting of eigenvalues 
of a reduced density p majorizes the probability distribution 
p(a) defined by the other reduced density a. In particular, 
the reduced density p strongly majorizes another reduced 
density a, written p >- a, if p(p) y p(a) and if the eigen- 
vector corresponding to p{p)\ coincides with the eigenvector 
corresponding to p(c)j. That is, this condition requires that 
there exists a common basis diagonalizing p and a. For 
more information about majorization, please see Bhatia's text 
book[16]. Using these notations, we can describe Nielsen's 
condition for LOCC transformation as follows. 

Lemma 4: Nielsen[4] We can transform an entangled state 
$ to another entangled state by LOCC if and only if a y p, 
where p (a) is the reduced density (partially traced state) of 
$ (40, respectively. 

Therefore, by using the above Nielsen's Lemma, the optimal 
performance of DFLEC, i.e., the maximum fidelity can be 
evaluated based on majorization as follows. 

Lemma 5: Let a be the reduced density of a given pure 
state \&, and p be the reduced density of the given initial pure 
state $. Then, we have 

max(*|C($)|*) =max max (tt y^^uY , (30) 
c p'hp ^unitary V > 

where the quantum operation C runs over all quantum opera- 
tions with LOCC in the maximum of LHS. If 4 is a maximally 
entangled state with the size L, i.e., the operator T :— La is 
a projection with the rank L, then the relation 

fJ7T) 2 



max max 
p'hp t/:unitary 



Tr VpVctlT 



Tr 



max ■ 

p'hp 



(31) 



holds. 



Proof: For any pure state VP, $, we have 

where two unitaries U\ and U2 are defined as 

U lP U* = Tm A U 2 aU; = Tt Ha 

Using Lemma|4] we can prove d30i . Next, we choose normal- 
ized basis {ei}f =1 and {fi}f—i as 

L 

i=l 

Using Schwartz inequality twice, we have 

L 

Tr^TU = Y i {fi\^/p i \^) 



\ i=l 



Since 

L L 

Y,{MV?\fi),Y,{ei\V?\ei) < max TiV^V*T, 
~[ ~[ y : umtary 

we obtain 

max TtV^/p'V*TU= max TrV^V*T. 
(7,V:unitary Viunitary 

Therefore, the equation 

max max Tr \[~p'TU = maxTr \[~p'T (32) 
p'hp t/iunitary p'hp 

holds because U pU* > p. Equations J30l > and d32l > guarantee 

CD- ■ 

However, it is not easy to directly connect the above lemma 
to the information spectrum. Hence, we prepare the following 
lemma for the evaluation of the RHS of ( I31> . This lemma plays 
an important role in the converse part of the main theorem. 

Lemma 6: When a projection T and an integer M satisfy 
TrT > M, and the two reduced densities p' and p satisfy 
p' y p, the inequality 

Tr y/^T 



(33) 




holds. 

Proof: Assume that TrT = N{> M). Without loss of 
generality, we can assume that p' y p. Let us diagonalize p 
and p' as p = J2i s i\ e i) ( e i\ and P' = J2i s i\fi)(fi\^ where 
Si > s i+ i,s- > s' i+1 . The inequality Tr y/p'T < Y,iLi V^i 
holds. We define the probability distribution {si^} and ijy as 

f N 1 

{si, N } := argmax I V ^ 



1 

N 



> s 



ijv + l ■ 



Similarly to i^, we can define %m- Since the function x 1— > y/x 
is concave, we can prove that S{ = s^jv for i < i^. Since 

im < In, 




<V«m 



i,^ g « = ^Tr |p > ^)^p{p 

N 

^/sZn < \/N -i 



\ i=iM+l 



= \/A r 



I'M 



i 



i=l 




Thus, we obtain (I33i . ■ 
In order to treat PFLEC, we have to consider a measuring 
operation with LOCC. Lo and Popescu characterize a pro- 
jection valued measure (Every P u is a projection, and 
Pu> is the identity.) on the system B as follows. 
Lemma 7: Lo and Popescu[17] For any projection valued 
measure {Pu,b} on the system B, there exist a projection 
valued measure {Pu>,a} on the system A and local unitaries 
U Ut A and Uu,b sucn that 

(/®P u ,b)|*) = (C^.ASt^.fl)^^® (34) 
That is, if the initial pure state is known, the operation 
corresponding to any projection valued measurement on B 
can be replaced by a projection valued measurement on A 
and local unitaries on A and B based on measuring data. 
However, we have to treat a general measuring operation with 
LOCC. The above Lo and Popescu's result can be generalized 
as follows. 

Lemma 8: Given a measuring operation I — {/^} with 
LOCC on a tensor product space Ha <8> Hb and a pure state 
|$)($| on the tensor product space Ha ® Hb, there exist a 
POVM {M w } (Every M w is a positive operator, and ^ w M w 
is the identity.) and the quantum operation C u with LOCC, 
such that 

I u ($)=C u (y/li^®I\$){§\y/Mu®I), Vcu. (35) 
Proof: It is known that any measuring operation Is = 
{Iui,b} on the system B can be described by the projection 
valued measure {P Uj _b} on an extended space H' B D Hb and 
quantum operations C^,b on B such that 

I<*,b(p) = Cu,b(Pu,brPu,b)- 
Applying 1341 . we have 

{Iu,B ®/)($) 

)((^,A®C/^)(^,A®/)|$)<$| 

(P w ,A®/)(^,A®^,b)*)- 

Hence, any operation on P can be described by the combina- 
tion of the projection measurement {P^^a}^ on A and local 
operations based only on the measuring data of {P u ,a]ui- 

Now, we focus on a measurement operation I' = on 
a tensor product space Ha <8> consisting of LOCC and a 
pure state |$)($| on Ha ®Hb satisfying the condition (A): 
the set fl = {uj} consists of all sent classical informations. 

Then, there exist the projection valued measure {Pu.b}u on 
an extended space H' B D Hb and quantum operations C w a 
and C w .s such that 

- (c w , A ®c UtB )({p UtA ®i)\®){®\{p^ A ®i) 

Even if the measurement operation / = {Ik]k with LOCC 
does not satisfies the condition (A), there exist a measurement 
LOCC operation /' = {I^^en with subset C O satisfying 
the condition (A) such that 



eigen value 



h 



(36) 



h(x) 



O 

Fig. 1. Illustration oi h(x) 

Hence, we have 

= E ic u ,A®c u , B )((p ll> , A ®i)\*)(9\(p u , A ®r) 



That is, there exist a projection valued measure {Pk, A } on an 
extended space H' A D Ha and LOCC operations such that 



Ik{$) = C k ((P k! A®I)\®)($\(P kl A®I) 



Since the projection Pu A to Ha satisfies that Pu A Pk,APi-L A — 
{Pk,APn A )* Pk,APn A > there exists a unitary U k , A such that 

KaPu a = \[m£ ■= Uk,APn A Pk,APn A - 
Hence, we obtain 



M^®/)|$)($|( A 



'M, A 0/)(C/ fe ,A 



J)' 



Therefore, the proof is completed. ■ 
In order to use the information spectrum method, one may 
characterize the optimal failure probability based on Tr p{p — 
x > 0} for the reduced density p of the initial state. However, 
it is difficult. Hence, we focus on h(x) := Tr(,o — x){p— x > 
0} instead of Tr p{p—x > 0}. Suppose that we wish to reduce 
all eigenvalues of the reduced density p to be no greater than 
x. This incurs a probability of failure given by h(x). Upon 
success we obtain a normalized state whose largest eigenvalue 
is not greater than x/(l — h(x)), which is majorized by a 
maximally entangled state of the dimension [(1 — h(x))/x\. 
It turns out that this method is optimal among PFLECs as 
follows. 

Lemma 9: The bound on the performance of PFLEC based 
on $ is evaluated by using the function h(x), as follows: 



max 

/={/ ,/ 1 }: PFLEC of 



{L(7)|Tr7 ($) < h(x)} 



-(l-h(x)) 



(37) 



where [^J denotes the maximum integer n satisfying n < x. 
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Proof: From Lemma [8] for any PFLEC I, there exist 
two quantum operations Co and C\ with LOCC and a positive 
operator P such that < P < I and 



Trl ($) 



Tr(J - P)p 
TvPp 



= c^ivT^ 5 ® r)\<p)(<$>\{vT^p ® i)) 

/„($) = Co((VPO/)|$)($|(\/P® /)). 

Hence, we obtain the following equations for the following 
reasons. 



linn <Tr/ (<l>) 
/={/ ,/i}:PFLEC of * 1 



Tr/ X ($) 



min {Trp(/-P)b- VPpV / P>0} (38) 
7>p>o on H 

/>P> T1 o i on « {Tr(p _ vWp)I* - VPPVP > 0} 09) 

min |1 — Tr clx — a > 0, p > a\ 
a on h ~~ ~~ 



(40) 

ei|cr|e !; ) < Si,x| (41) 

1— Si— x = Tr(p — x){p — x > 0} = h(x), 

(42) 



= min < 1 — > (e,'|cr| 
CT on w | ^ x 1 1 



i:Si >x 



where we diagonalize p as p = X^i s il e i)( e il m <i41i . From 
Lemma[4] there exists a quantum operation Ci with LOCC that 
transforms the state ^^-(P®7)|$)($|(P®/) to a maximally 
entangled state with the size L if and only if Tr J j.W > PpP. 
Thus, from Lemma |SJ we obtain 1381 . In general, for any 
bounded operator A, there exists a unitary operator U such 
that AA* = UA*AU*. Thus, the condition x - PpP > is 
equivalent with the condition x — ^JpP^fp > 0. We obtain 
( 1391 . Replacing ^fpP^fp by a, we obtain J40I . 
Equation (I42> implies 



max {L(/)|Tr/ (a>)<^)} 
/={/ ,/i} : PFLEC of * 



= max < —(1 — h(x)) 



-^■(1 — h(x')) is an integer, 

h{x') < h{x) 



~(l-h{x)) 



where the second equation follows from the fact that the func- 
tion h(x) strictly monotonically decreases and is continuous. 



VI. Asymptotic theory 

In this section, based on non-asymptotic formulas given in 
section [V] we prove our main theorem. For this purpose, we 
need to prepare the finite-version of the information-spectrum 
quantities for a projection operator T„ and a reduced density 



it f\t\ r—i . qc Tr\llr\vi7c 
(J n Ull /1 1 j± yi lUUUWs. 




C (T \a ) 

Sn V n 1 u n ) 


— _ A. loo-TriT T 
n 


Sn,l/2\ ± n\<Jn) 


n 


n (T \ 


— — — lncr Tr( T T \ 

— — log ir J l n , 

n 




= --logTra n {I-T n ), 
n 


Cn,l/2( T nWn) 


= -i log TV V5^(J-T n 



As the limiting version, we define 



C(T|<x) 


= ]im( n (T n \a n ), 


C(f\ff) 


= limCn(r„|cr„), 


Ci/ 2 (T|c?) 


= lim( n (T n |<7 n ), 


c 1/2 (?V) 


= limCn,i/2(r«|cr n ), 


Tj(T|*) 


= lim?7„(T„|cr„), 




= lim??„.(T»|g„.), 




= ISlC(r„k n ), 


C=(f|<x) 


= }m(n( T nWn), 


Ci/ 2 (T|<?) 


= limCn(r„|cr„)) 




= limC,i/ 2 ( r «|crn), 



for sequences cr = {a„} and T = {T n }. For the 
projection S n (a) := {p„ < e~ na }, we simplify 
Cn(5 , n(a)k„),Cn4/2(' S 'r l (a)|cr„), ?7„ (S„ (a) ) , C£ (^n (a) |ct„ ), 

and C ! i/ 2 (' S '™( a )l (7 n) t0 Cn(a\<T n ),Cn,l/2(a\(Jn), 

Vn(a),(n( a Wn), and Cli /2 ( a l cr ™)- We can similarly 
define ((a\v), ((a\a), C 1/2 (a\a), C x/2 {a\a), f}(a\a), v(a\a), 

/ 2 (a\<?), and C[ /2 (a\^)- Using these 
values, we can characterize the RHSs of J3 II . (I33i and 
(I37> . In particular, when a sequence <r equals the sequence 
P = {Pn} °f me reduced density of the given state, we omit 
in the above values. 

Moreover, to discuss the asymptotic theory, we need to 
define the concept "majorization" in regard to sequences of 
reduced densities. The sequence of reduced densities er = 
{<7„} majorizes (strongly majorizes) another one cr' = {cr' n }, 
written er >z cr' (cr >~ cr') if a n > a' n (a n >- a' n ), respectively. 

In the following, we proceed to the proof of our main 
theorem. Before it, we should remark that in an asymptotic 
case, we can neglect the gap between [£nj and L n because 
L n is large enough. 

Lemma 10: Without any assumption, the equations 

Bi(e) = B 2 {e) = sup{R\K(R) < e}. 

R 

hold for every e £ [0, 1]. 

Proof: From the definition, the inequality B>i(e) > B 2 (e) 
is trivial. We only need to prove the two inequalities 



B 2 (e) > sup{R\K(R) < e} 

R 

Bi(e) < sup{R\K(R) < e}. 

R 



(43) 
(44) 
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Let R be a real number satisfying 

K{R) < e. 



(45) 



From Lemma [9] there exists a PFLEC I n such that 

Tr = h n {e~ nR ) and L n = e nR (l-/i n (e- nR )), where 
h n (x) := Tr(p„ — a;){p n — a; > 0}. From ( 145 1 . we have 

lim - log L„ = R, 

n 

lim" Tr /"($„) <E5Trp n {p„ - e -" fl > 0} 
< e. 

We have now obtained the direct part (I43> . 

Next, we proceed to the converse part (1441 . Let 7™ be a 
DFLEC satisfying lim(*„|C" ($„)|*„) > 1-e. For any i? < 
lim — log L n , we have 

e nR 

lim — = 0. 

From Lemma |6] for any T n satisfying TrT„ = L n , we have 
(Tr V&) 2 



<— ( x/Tr{p„ > e-^V^Mpn > e ~ nii } 



in - Tr{p„ > e-™ R }JTr p n {p n < e~™ R } 



ITr{p n > e- nR } 



Trpn{Pn > e~ nR } 



1 Vl* p n {p n <e-» R } 



Since lim ^*i£n=S. 1 < li m 

,. (Tr^Tj 2 



0, 



1 — e < lim ■ 



< limTrp„{p„ < e- nH } 



=1 - K{R). 

Thus, we obtain (I44i . ■ 
Lemma 11: We have 

B e , D (r) = B e , P {r) = sup{R\{ c (R) > r}. 

R ~ 

Proof: Since B e ^o{r) > B B: p(r), we only need to prove 
the inequalities 



B e , D {r) < su P {i?|C c (i?) > r} 

R ~ 

B e . P (r) > su P {i?|C c (i?) > r}. 

R ~ 



(46) 
(47) 



First, we prove the direct part J47I . Assume that ( C (R) > r > 
0. From Lemma [9] for any R, there exists a PFLEC I" with 

the size e nR {l - (1 - i„(i?))e-<^)) such that 



Tr^(<f> n ) = (l-t„(i?))e-" C " (fl) , 



where 



t„CR) := 



-nR 



Tr{p n > e~ nR } 
Trp n {pn > e~ nR } 



Since C c (-R) > 0, we have < (1 - t n (R)) e - n ^ R ^ < 
e ~ n Q(R) o. Thus, we have the following relations 

limiloge" fl (l- (l-i„(i?))e-™ c " (fl) )) = R 

lim — log Tr > ( C (R) > r, 

n — 

which imply the inequality J47I . 

Next, we proceed to the converse part ( 1461 . Assume that the 

DFLEC (C n ,*„) satisfies 



lim- log (1 - (¥„ | C n ($„)|*„)) > r. 



(48) 



We define the projection T n and the reduced density as 

(Tr J~$T n 

T n := L n Tr HB |* n )(*n|, /»„ := argmajy^ 
Then, Lemma [5] and ( I48> yields that 



L n 



lim — log 1 
n \ 



(T rv ^r„) 2 



> r. 



For any R' < R := lim — logL„, there exists an integer 
N such that R n := ^logL™ > R' for Vn > AT. When a 
projection T„ satisfies that TrT„ = L n , Lemma [6] implies 
that 

(Tr^^) 2 
L n 



< —[ \/ Tl {Pn > e nR '}\/Trp n {p n > e nR '} 



L n - Tv{p n > e nR '}^Trp n {pn < e nR '}y 

<^ e -?('7n(H')+C(R')+-Rn) 

+ y/l - e -n(«„(fl')+-R„) v /l _ e -"CS(fl')^ 2 
£^1-1 (e-n(.T,n(.R')+Rn) + e -«C(fl')j 

+ e -f (r h JR')+CAR')+R^y 



e -f )+«„) 



(49) 



Since e"t < e -f («» -•#)<;-§ ("»(«')+«') < 
e -f (fl ra -fl') e -f CSC-R') < e-fCSC-R'), we have 



^ e -fC«(«') _ e -f (77„(fl')+«n) 

>(1 - e -f (fl..-«'))2 e -«C(fl'). 



Thus, 



5 (■-*'* 



(R')+Rn) _ g-fC(ii') 



< fl - 1(1 - e~9 (R n - R ')f e -nQ(R'l 



(50) 



Since lim(l — e R ^) 2 = 1, it follows from d49l and 

® that 

C c (i?')> lim - log fl - (^V^ T n) 2 \ ^ r 
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Since R 1 is an arbitrary real number satisfying R' < Rq, the 
relation Rq < sup R {R\ ( C (R) > r} holds. Therefore, we 
obtain TO , _ ■ 

Lemma 12: When £(a) — £(a) —: ((a) and there exists a 
real number a such that ((a) < C c ( a )> 

= sup{a — min{C(a), a + min{C(a), a + fj(a)} < r} 

a 

= inf{a — min{((a), a + min{£(a), a + fj(a)} > r} 

a 

= inf{a-C(a)|C(a)<r}. 

a 

Proof: First, we prove the direct part. Consider a PFLEC 
J" satisfying 



in — 



l-/i n (e- no ) 



Tr/ n (4>„) = M e -" a ). 



Thus, we have 



lim — logL„ =0 — min{£(a), a + r](a)}, 
n 



-1 



lim log (Tr /"($„)) = min{C(a), a + 77(a)}. 

n 

Therefore, we have 

> sup{a — min{£(a), a + /7(a)}| min{£(a), a + rj(a)} < r} 

a 

= max < sup{a — £(a)|£(a) ^ r }i snp{—rj(a)\a + 77(a) < r} 

^ a a 

= sup{a- C(a)|C(a) < r}, 

a 

where the final equation is derived by Lemma IT31 as follows. 
Using Lemma [TBI we have sup a {a — £(a)|£(a) < r } > 

su Pa{-!Z( a )l a + !Z( a ) - 7 '} - SU Pa{~V( a )\ a + V( a ) < 
Next, we proceed to the converse part. Let {/"} be a 

sequence of PFLECs such that r > lim — log(l — e„), where 

e„ := Tt/q (<!>„). In the following, we focus on lim — logL„. 

Let a be a real number satisfying 

lim — - log(l — e„) < r < min{C(a), a + 77(a)}. (51) 
n 

Since 

DS — log (Tr p„{p„ < e-" a } + e-" Q Tr{p„ > e"" }) 
71 

= min{C(a),a + 77(a)}, 

there exists an integer N such that 

Trp„{p„ < e~ na } + e~ na Tr{p n > e~ na } 
>l-K{e" na ), Vn>N. 

Lemma |9l guarantees that 

e na (Trp n { Pn < e- na } + e - na Tr{p n > e -" a }) 



1 - h„ (e 



> L„ 



(52) 



Taking the limit of the exponent, we have 

lim — logi n < a — min{£(a), a + 77(a)}. 
n 



From i51\ . we have 

Bt P (r) 

< inf{a — min{C(a), a + 77(a)}! min{£(a), a + 77(a)} > r}. 

a 

It follows from (I70i that the function a 1— > min{C(a), a+7j(a)} 
is continuous. Thus, 

inf{a — min{£(a), a + 77(a)}! min{£(a), a + ?7(a)} ^ r } 

a 

= sup{a — min{£(a), a + 77(a)} | min{£(a), a + 77(a)} < r}. 

a 

The proof is now completed. ■ 

Lemma 13: When £(a) = C( a ) — : C( a ) an d there exists a 
real number a such that £(a) < C c ( a X 

= sup sup{-limry„(r„)|lim2C„. 1/2 (r„|p n ) - 7j(T n ) < r} 
p'hp f 

(53) 

= sup{-limry n (T„)|IirK2^ )1/2 (T„) - ry(T„) < r} (54) 



T 

= sup \a — r 

sup <; I - inf I ((a') - |- 



inf C(a') " 



a! < a \ + — < r 



(55) 



a' < a 



inf CK) 



a' < a K g < r J> . (56) 



Proof: Equation J53) follows from yTJ. Since the 



1' < a i 



is continuous 



function a 1— > inf a / < C( a ') — % 
and decreases monotonically and the function a 
inf a ' I C( a ') ~ t a ' — a | +§ is continuous and increases 
monotonically, equation d56l > holds. First, we prove the direct 



part: 



sup{-limr? n (T n )|lim2C 1/2 (T n )-r7(T n ) < r} 



> sup ^-inf{ C (aO-| 
inf (((a') - 



a' < a 



a' < aj + - < rj . (57) 

As we prove later, we can choose a projection T n (a, R) such 
that 

r] n (T n (a,R)) = -R, (58) 
Q tl/2 (T n (a,R)) < max{C il/2 (a),- J R+ |} . (59) 

When f] n (a) > -R, the projection T n (a, R) := {p n -e~ na > 
0} satisfies ( I59> . Otherwise, the projection T n (a, R) := {p n — 
e -na > o} + _ e -na < } - f n (a,R)) satisfies ^J, 
where T n (a,R) is constructed as follows: We choose m := 
e nfl normalized eigenvectors {e^}™ j of {p„ — e~ na < 0}p n 
in descending order concerning the eigenvalue, and define the 

projection f n (a,R) by YZi \ e i)«\- The choice of { e ii^i 
and the relation e nR = Tr{p n - e - na < Q} e - n( -- R - r i" ( - a '>') 
guarantees 

TrVP^iPn ~ e- na < 0}e- n (- R -^ a » < Tr ^f n (a,R). 

(60) 
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Then, we can check the condition ( I59> as follows: 

Vn {T n {a,R)) = —\ogTr(I-T n {a,R)) 
n 

=—logTif n (a,R) = —lo ge nR = -R, 

n n 

C,i/ 2 (T n {a,R)) = — log Tr v /p^T„(a, R) 



1 



-n{-R-r, n (a)) 



<_lo g TrV^{p„-e- na <0}e 

n 

=Cn,l/M -R-Vn{a) <-R+\- 

Now, we apply Lemma [n] to the case p n — p n , cr„ = y/p^T- 
Since 



{p„ - e na a n > 0} = {p r , 



P« > 0} 



e 

={VP^ - e" a > 0} = {p„ - e 2 " a > 0}, 

we have 

V(a) = (_1 /2 (MP), C(a)=C(2a|p) 
From Lemma El the maximum a r of 

a 



(61) 



a' < a ^ + 



± inf C(a') - - 

2 a' 1 2 



a' < a T 



inf|C(«') 
exists. We define R by 

i? 

Then, i? equals to the right hand side of ( I57K and we have 

K^2Q A/2 (T n {a k ,R)) - r)(T n {a k ,R)) 

<2 max |d/ 2 ( a i- + V&)> 



1/fc 



where a k ■= a r + 1/fc and fc is a fixed integer, and the last 
inequality follows from J84l > in LemmaEl m Appendix|X| We 
define N k as the minimum integer satisfying 

2Qi/2( T n{ak, R)) - r){T n {a k ,R)) <r+^, Vn > N k . 

For the sequence b n :— mm k {ak\n > N k }, we have 

^2Q tl/2 (T n (b n ,R)) - ri{T n (b n ,R)) < r. (62) 

Inequality d57i follows from (I62> and the first equation of ( 1591 . 

Next, we prove the converse part. Assume that {(T n ,p' n )} 
satisfies lim n _ l . 0o 2^ , 2 (T n \p' n ) - r)(T n ) < r. There exists 
a subsequence {n k } such that lim r\ nk {T Uk ) = — Rq := 
lim 7?„(r„). Focusing on the projection {p' n — e~ na > 0} = 
{\fPn~ e- na ' 2 > 0}, we have 

Tr ^K{Pn ~ e~ na > 0} - e na ' 2 Tr{p' n - e'™ > 0} 
> Tr v^(J - T„) - e" Q / 2 Tr(7 - T„), 

which implies 



>TrV^(/-r„). 



>0} 



na/2 



Tr(/ - T n ) 



Taking the limit fc — > oo, we have 

min ) J im Q.i/aOKJ' | - }• < lim C fc ,i/ 2 ( T nJP«J- 

[ n, f OO ^ j ft- t OO 



Now, we apply Lemma IT4l to the case p Tl 
this case, similarly to 16 H . we have 

2(«) = C C 1/2 (2a|p0. C(o) = C(2a|p')- 
Hence, {73) yields that 

J^c, 1/2 («ip;j>c; /2 («ipo 



>inf <j C(a'\p'^ 



a' < a 



Since p' n y p n , we have ((a'\p') > ((a 1 ) = ((a r ), i.e., 



MU(a'\p')-- 



'<a\>m£\«a')-~ 



a < a> . 



Thus, 

r> lim" 2Q !l/2 (T n \p' n ) - r)(T n ) 

n — >oo ' ' 

> ns 2Q >1/2 (T„ fc | P ;j-7j(T„j 



2llml ^ fe !™,^fc.l/2( a l^nJ> 2 ~ Ra 

>2min <{ inf i ((a 1 ) - 



a' < a 



a 
2 



Ro 



Since the function a i— > f — inf a ' | £(a') — %■ 
continuous, there exists a real number a such that 



- Ro- 
(63) 

a' < a j is 



flo = f-inf C(a')-| 



a' < a 



Using d63i . we have 

r > inf <j C(a') 

a' 

which implies 



i?o<sup<j--inf<j C(a') 



a'<aj. + -, 



a' < a 



inf<j C(a')- 



a' < a ? + 2 < r 



The proof is now completed. 



VII. Relation to random number generation 

As a related problem, it is known to transform from a 
given known probability distribution p to a desired probability 
distribution q. If it is possible, the majorization relation q >z p 
holds. However, even if the majorization relation q y p holds, 
this transformation is not necessarily available. Hence, if the 
two entangled pure states $i and $ 2 have Schmidt coefficients 
corresponding to p and q, the Quantum LOCC operation 
transforming from $! to $ 2 is easier than transform from 
p to q. 

In particular, when the desired distribution is the uniform 
distribution, this problem is called intrinsic randomness. In this 
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problem, our operation of intrinsic randomness is described by 
the map ip from the original space ft to M. = {1, . . . , M}. 
When the initial distribution is p and the uniform distribution 
is described by puj on A4, one of criteria of its quality is the 
half of the square of Hellinger distance between pot/j^ 1 and 
Pm' 



M 



e(ijj,p) :=l-^2> 



M ' 



(64) 



In this case, we describe the size of its target uniform dis- 
tribution ip by M(tp). Hence, for a sequence of the initial 
distributions {p n }, we can define the optimal rates 



BeAr) 

b;. h {t) 



sup < lim 
sup < lim 
sup < lim 



logM(^„) 



n 

log M(ip n ) 



n 

log M(lp n ) 



\ime(ip n ,p n ) < ej 

Urn — loge(^ n ,p„) > r\ 
n ) 



n 



lim — log (1 - e(ip n ,p n )) < r\. 
n J 

The variational distance version with the constant constraint 
has been discussed by Vembu & Verdu [19] and Han [13]. 

Let $„ be the entangled pure state with the Schmidt 
coefficient corresponding to p n . When C„ is the quantum 
LOCC operation corresponding to ip n and ^> n is the maximally 
entangled state with the size M{ip n ), we have 



1 



:(lp n ,Pn) = V(*n|Cn( $ n)|*n), 



(65) 



As is shown Hayashi[22], the relation 

Bkl(£) = sup{a - C(a)|C(a) < e l 

a 

holds. When ("(a) is continuous, 

B KL (e) =B* eiP (e). (67) 
In particular, if the limit of Renyi entropy is differentiable, 

B KL (e)>B* eH (e/2) (68) 

when e < — hip (|) — ip (¥)■ The above relation is an inter- 
esting relation between Hellinger criterion and KL divergence 
criterion. 

VIII. Concluding remarks 

We derive asymptotic bounds based on several formulations 
from Lemma [5] |6] and [9] Since these bounds are tight in 
a general source, the evaluations given in Lemma |5j |6j and 
|9] are useful in a non asymptotic case as well as in an 
asymptotic case. Even if the class of DFLEC is wider than 
that of PFLEC, their asymptotic performances are almost 
equivalent. A difference appears only between B* D (r) and 
B* P {r). For example, when the limit of Renyi entropy ip(s) 
is differentiable, B* D (r) is larger than B* P (r) if and only if 

r is greater than — hip (h) — ip (h). From d54l of Lemma [T3l 
the bound B* D (r) can be attained without an LOCC, i.e., the 
original reduced density p n is close enough to an appropriate 
MES only in regard to B* D (r). As a byproduct, in Appendix 
1X1 we establish several general relations between information- 
spectrum quantities. 



2e(ip„,p n ) - e{ip n ,p n f = 1 - (iqC n ($„)|tf n ). (66) 

Hence, comparing the entanglement concentration with the 
initial entangled state $„ and the intrinsic randomness with 
the initial distribution p n , (I66> yields that 



Since 



B H (e) <B D (2e-e 2 ). 



e(ipn,Pn) < 1 - (*„ | C„ ($„)|*„) < 2e(iP 

n : Pn ) : 

the inequality 

BeAr) < BeAr) 
holds. Moreover, the equation d65i yields that 

B* e A r ) ^ B l D (2r). 
When we adopt the KL divergence criterion: 

D(p M \\p°ip- X ) :=logA/ + ^^log ( P«J» 

i=l \wGi/'- 1 (i) / 



we focus on the following value: 

logM(^„) 



Appendix 

A. General relations for information spectrums 

Here, we prove some lemmas required by our proof. In 
this section, we treat information-spectrum quantities with 
more general definitions, which are given in Nagaoka and 
Hayashi[14]. This is because we need such a general treatment 
in our proof of Lemma [T3l 

For the two sequences {p n } and {cr„} of trace class positive 
semidefinite operators, we discuss how to characterize an 
information-spectrum quantity r/(a) := lim — log Tr a n {p n — 
e ™ a Cn > 0} by using two other information-spectrum 
quantities ((a) := lim^- logTrp,^^ - c na (j n < 



0} and ( c (a) := lim — logTrp„{p r , 



> 0}. 



As discussed later, when £(a) := lim — log Tr p n {p n — 
e~ na a n < 0} equals ((a) for any a, we can use the same 
method to characterize another information spectrum 77(a) := 
lim ^ logTrcr n {yO n — e~ na a n > 0}. As was proven by 
Nagaoka and Hayashi[14], the function £(a) increases mono- 
tonically, and other functions C c (°0 an d 77(a) decrease mono- 
tonically [14]. Focusing on the projection {p n — e~ na a n > 0}, 
we have 



v„ > 0} > 0, 



Bkl{<P) ■= sup <^ lim 

{^n} L 



Tr(p„ - e" 
which yields to 

Hm£)(p M(V)n) ||po ip- 1 ) < e} Tlpn { pn _ e ~ na a n > 0} > e - ,la Tr a n {p n - e-"V„ > 0}. 
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Thus, we have 

C C («) <v(a)+a. (69) 

Similarly, we can prove 

Tr{p n - e-™a n ){p n - e- na a n > 0} 
> Tr( Pn - e- na a n ){p n - e~ nh a n > 0}. 

By adding e~ na Tra n to both sides, we have 

Trp n { Pn - e- na a n > 0} + e~ na Tta n {p n - e - Ba oy, < 0} 
> Tr p n {p n - e- nb a n > 0} + er na Tr a n {p n - e~ nb a n < 0}. 

Taking the limit n — ► oo, we obtain 

min{C(a), a + 77(a)} > min{C(o), a + v(b)} (70) 

for any a and o[14]. When £(a) = £(a) for any a, we can 
replace 77 by 77. From inequality (170b . We can derive the 
following two formulas; 



77(a) + a > ((b) if 77(6) > 77(a) 
C(a)>a + ?7 (6) ifC(a) < C(b), 



(71) 
(72) 



which play important roles in the following lemmas. As a 
lower bound of 77(a), the following lemma holds. 

Lemma 14: If there exists a real number ao such that 
C(°o) < C C (°o!h the relations 



77(a) > inf{C(a') - a'\a' < a} 

— a' — 

= inf{C(a') -a'\a' < a} 



(73) 
(74) 



hold. 



Proof: From (1691 . the relations 

((ao) < C c (ao) < ?7_(ao) + ao 
hold. Since 77(00) > C( a o) — a o> we nave 

??( a o) ^ m f{C( a ') ~ a'|a' < ao}. 

— a' — 

For any a < ao, the relation ((a) < ( c (a) holds. Since ((a — 
0) < ((a), the equation (I74t holds. Similarly, we can prove 
that a real number a(< ao) satisfies (1731 . 

Next, we prove d73l for any a > ao by the transfinite 
induction. Assume that the relation d73l holds for any real 
number b satisfying a > b and 

77(a) < inf{C(a') - a'\a' < a}. (75) 

— a' — 

For any e > 0, we have 

77(a) < inf{C(a') — a'|a' < a — e} < i](a — e). 

— a' — — 

From 17 H . we have 77(a) > C(a — e) — a. Since e is arbitrary, 
we obtain the inequality 

77(a) > inf{C(a') - a'\a' < a}, 

— a' — 

which contradicts assumption d75t . ■ 
The following lemma is another characterization of the 
lower bounds of 77(a). 



Lemma 15: We obtain the inequality 

sup{a — C(a)|C(a) < 7'} > sup{— 77(a) |a + 77(a) < r}, 

a a 

which is equivalent to another inequality 

inf{C(a) - a\((a) < r} < inf {77(a) |a + 77(a) < r}. 

a — — a — — 

Proof: We prove it by reduction to absurdity. Assume 
that there exists a real number ao such that 

a + 77_(a ) < r, (76) 
-77_(a ) > sup{a - C(a)|C(a) < r}. (77) 

a 

We will lead contradiction with the two cases, case 1: a\ := 
inf a {a|?7(a) = 77(00)} > ao, case 2: ai = ao- 

In case 1, for any real number e £ (0, 00— ai), the inequality 
77(01 — e) > 77(01 + e) holds. Using d7 It . we have 

C(ai — e) < 77(01 + e) + ai + e = 77(00) + ai + e 
<r + (ai — a ) + e < r. 



Thus, 

sup{a — C(a)|C(a) < r} > ai — e — £(ai — e) 
■ n 1 — e — (ai + e) - 77(01 + e) = — 77(00) — 2e. 

Taking the limit e — > 0, we obtain supja — C(a)|C(a) < r} > 
—77(00), which contradicts d77i . 

In case 2, the inequality 77(00) < 77(00—6) holds for Ve > 0. 
Using i71\ , we have £(ao — e) < 77(00) + a < r. Thus, 

sup{a - C(a)|C( a ) < t-} > a - e - ((a Q - e) 

a 

>a - e - a - 77(a ) = -e - 77(00). 

This also contradicts d77i . ■ 
Define the sets / and /' as 

/ := {o e R|C(o) > ((a - e) Ve > 0}, 
I' := {a e R|C(a + e) > £(«») Ve > 0}. 
As upper bounds of 77(a), we have the following two lemmas. 



Lemma 16: We have two inequalities 

77(a) < inf {((a ) — a' I a' < a}, 
77(a) < inf {((a') -a'|a' < a}. 



(78) 
(79) 



If C(a) = C( a ) f° r an y rea l a ' we nave two other inequalities 



77(a) < inf {((a') - a' I a' < a}, 

a£l — 



(81) 



??( a ) < m f {C( a ') — a'|a' < a}. 
Proof: First, we prove d78b . Let a' G / be a real number 
satisfying a' < a. From (1721 . we have 

o'-e + 77(o') <C(a'-e), Ve > 0. 

Since e > is arbitrary, we obtain the relation 

77(a) < 2 (o') < C(o' - 0) - a' < C(a') - a'. 

From the arbitrariness of a', the above relation implies ( I78l l. 



Similarly, we can prove d80b . 
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Next, we prove ( I79> . Let a' G I' be a real number satisfying B. Gartner-Ellis theorem 
a 1 < a. From (72), we have Here, for our proof of TheoremE] we discuss Gartner-Ellis 

theorem [20]. Let X n be a sequence of random variables. 
Then, the logarithmic moment function is defined as 
If e > is small enough, K{f) . = ^^y^ 



a' + rj{a + e) < C( a ')- 



"(a) < v( a ' + e ) < C( a ') — a ' ■ 



where Ex„ denotes the expectation concerning the random 



_ , . . . , , , . ..... jssjTi variable X n . The logarithmic moment function A n (t) is con 

From the arbitrariness of a , the above inequality implies i/9t . 



Similarly, we can prove d8 1 1 

Lemma 17: Assume that a real number r satisfies that 



vex. 



Theorem 18: Assume that the limit A(t) := liirin^oo 
exists. Then, defining the rate function 

r <sup{M{((a')-a'\a' <a}+a}. (82) A* ( R) := sup tR- A(t\ (88) 

a t 
The maximum a r of we have 

{a\mf{C(a')-a'\a'<a} + a = r} (83) 

exists. Moreover, the inequality 

?7(a r + e) < inf (C(o) - a\ a < a r } , Ve > (84) 

holds. When £(a) = ((a) for any a, we can replace 77 by rj in 
the above argument. 

Proof: Since the function q : a 1— ► TT . , , , , - , ,. ,, . 

. r / ' ,1 , -, . . 3 , . Using the above theorem, we can show the following theorem. 

mr„/ ( b -a a <a + a is continuous and increases „. , _ . . /jA . , . . 

— ■ ,, • ' r. ,, — ,- , . , , , Since the function A„(r) is convex, the A(t) is convex, too. 



lim 


-1 
n 


log 


p.] 


\x n 

I n 






inf A*(fl) 

i£>a 


(89) 


lim 


-1 
n 


log 




\x n 

I n 






inf A*(R) 

R>a 


(90) 


lim 


-1 
n 


log 


p.] 


[x n 

I « 






inf A*(iZ) 

it<a 


(91) 


lim 


-1 
n 


log 


p^ 


fA„ 
I n 


< aj 




inf A*(R). 

R<a 


(92) 



Hence, when we choose the real numbers R\, i?2, -R3 and i?4 



as 



not belong to the set (I83L the relations 

£(a) < C( a r) = inf { C( a ') — a'| a' < ar} + <^r = ^ 
< inf { C(a') - a \ a' < a r + e} + a r + e < ({a r + e) 



monotonically, it follows from (I82> that set (I83i is bounded 
and closed. Thus the maximum of the set (1831 exists. 

Next, we prove (I84> . First we assume that ^ ^ 

C(a r ) - a r > inf { C(a') - a'| a' < a r } , (85) fll := t^™ ~T' ^ 2 := t^+o (93) 

n •= lim M o . = i im AM f 94^ 
Since the function g increases monotonically and a r + e does 3 ' t^_o i ' t ' 

the relations 

R4<Rs<R2< Ri (95) 

hold. Thus, as is proven latter, the equations 

hold for a < a r . Applying i72l to the case b = a r + e, we ,. -1, f X n 1 .. -1, _ f X n 

. __ rr ^ — lim — log Fx < > a > = lim — log Fx { > a 

obtain 184). n 6 "\n-J n 6 "\n 

Second, we assume the opposite inequality f q if a < R2 

\ . ■ ( r» (1 /i / ^ 1 =< maxti?- A(t) > if i? 2 < a < R\ (95) 

C(a, ) — a r < mi { C( a ) — a a < a r ) . (80) j t>o v ; 

_ a' - |^ cxd if i?i < a 

There exists a sequence {a„} such that anc j 

C(a„) - a„ - inf { C(a') - a'\ a' < a r } ^ -1 lQgp < a \ = lim ll l og ? Xn / ^ 

[ if R 3 < a 

From the above relations, there exists an integer N such that = / maxti? — A(t) > if R4 < a < R3 (97) 
C(a n ) < ((a r ), Vn > N. Using (7JJ, we have I oo^ if a < R 4 

rj_(a r ) < C( a »i) — a n- hold. Moreover, if the function A is differentiable at t > 0, 

and if R2 < a < A' (to), we have 



< a 



Thus, we obtain 

77 (a r ) < inf (C(a') - a'\ a' < a r } , (87) 

_ = sup tR-A(t), (98) 

which implies i84l . ■ t >t>o 



-1 \X n 1 -1 [X n 

lim — log P x \ >a> = lim — log Px s > a 

n In n In 
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Proof of and First, we calculated the rate 

function A* (a). When R 3 < a < R 2 , 

A* (a) = sup to - A(t) = Oa - A(0) = 0. 
t 

Assume that R% < a < R\. Then, if e > is sufficiently 
small, 

A* (a) = sup to - A(t) > ea - A(e) = (a - i? 2 )e + R 2 e - A(e) 
t 



(a- R 2 )e + R 2 e 



Mm M e 



(a~R 2 )e > 0. 



Now, we choose to 7^ such that t a a = A(t a ). The convexity 
of A guarantees that 

A* (a) = sup to — A(t) = max ta — A(t). 

t o<t<t a 

For a such that R 2 < a' < a, since t a > t a >, we have 

A*(a') = max to' - A(t) = max to' - A(t). 

o<t<t a , o<t<t a 

Hence, the function A* is continuous [R 2 , a] . Thus, the func- 
tion A* is continuous [R 2 ,Ri). in addition, when a > R\, 
A* (a) = 00. Hence, when a < R\, we obtain 



inf A*(R) = inf A*(R) 

R>a R>a 



max tR - 



A(t) > 



if a < R 2 

if R 2 < a < Ri 



When a > 



inf A*{R) = inf A*(R) = 00. 

Therefore, we obtain ( I96> . Similarly, we can prove ( 197k 

Moreover, for a such that R 2 < a < Ri, we choose i' Q = 
argmax t ta — A(t). The convexity of A guarantees that when 
R 2 < a' < a, we have t' a , < t' a . Therefore, we prove ( 1981 . ■ 
Finally, in order prove dl6l and (I17> in our proof of Theorem 
13 we focus on the probability distributions p n = {p n ,i}, and 
apply the above discussion to the random variable — \ogp n .i. 
Using J9fjj, (g7j and (|98}, we obtain (|l6j and O- 
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